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Abstract. The nature of the mechanism responsible for producing the spectacular, geometrically thin, spherical 
shells found around some carbon stars has been an enigma for some time. Based on extensive radiative transfer 
modelling of both CO line emission and dust continuum radiation for all objects with known detached molecular 
shells, we present compelling evidence that these shells show clear signs of interaction with a surrounding medium. 
The derived masses of the shells increase with radial distance from the central star while their velocities decrease. 
A simple model for interacting winds indicate that the mass-loss rate producing the faster moving wind has to 
be almost two orders of magnitudes higher (~ 10“® Mq yr“^) than the slower AGB wind (a few 10“^ Mq yr“^) 
preceding this violent event. At the same time, the present-day mass-loss rates are very low indicating that the 
epoch of high mass-loss rate was relatively short, on the order of a few hundred years. This, together with the 
number of sources exhibiting this phenomenon, suggests a connection with He-shell flashes (thermal pulses). We 
report the detection of a detached molecular shell around the carbon star DR Ser, as revealed from new single-dish 
CO (sub-)millimetre line observations. The properties of the shell are similar to those characterising the young 
shell around U Cam. 
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1. Introduction 

The late evolutionary stages of low- to intermediate-mass 
stars, as they ascend the asymptotic giant branch (AGB), 
are characterized by the ejection of gas and dust via a slow 
(5 — 30kms“^) stellar wind. This mass loss (on the order 
of 1O“^-1 O“®M 0 yr“^) is crucial to the evolution of the 
star and a key process for enriching and replenishing the 
interstellar medium (e.g., Schroder et al. 1999; Schroder 
& Sedlmayr 2001). Yet, the mechanisms by which it oc¬ 
curs are not understood in detail. In particular, there is 
evidence that strong variations in the mass-loss rate, per¬ 
haps related to flashes of helium shell burning (Olofsson 
et al. 1990), can lead to the formation of circumstellar de¬ 
tached shells. The existence of such shells can be inferred 
from their excess emission at far-infrared wavelengths (due 
to a lack of hot dust close to the star, Willems & de Jong 
1988; Zijlstra et al. 1992), and are confirmed by a ring¬ 
like morphology in maps of dust or CO line emission (e.g.. 
Waters et al. 1994; Olofsson et al. 1996; Lindqvist et al. 
1999; Olofsson et al. 2000). 
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Detached CO shells have been observed around about 
a half-dozen AGB stars (see review by Wallerstein & 
Knapp 1998), all of them carbon stars. However, only 
U Cam (Lindqvist et al. 1996, 1999), TT Cyg (Olofsson 
et al. 2000), and S Set (Olofsson et al., in prep.) have been 
mapped in detail with millimetre-wave interferometers. 
These three cases reveal remarkably symmetric shells ex¬ 
panding away from the stars. However, there are also clear 
indications of a highly clumped medium. High-resolution 
molecular line observations offer the possibility not only to 
confirm detachment, but also to study variations in mass 
loss on short (<100 yr) timescales, related to the shell 
thickness, and to look for departures from spherical sym¬ 
metry. Single-dish CO maps provide evidence for detached 
CO shells also in the cases of R Scl, U Ant, and V644 Sco 
(Olofsson et al. 1996). 

Gonzalez Delgado et al. (2001, 2003) have imaged the 
circumstellar media of R Scl and U Ant in circumstellar 
scattered stellar light. They found evidence for both dust- 
and gas-scattered light in detached shells of sizes compa¬ 
rable to those found in the CO line emission. These data 
provide high-angular-resolution information, but the sep¬ 
aration of dust- and gas-scattered light is problematic. 
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Based on colour-colour diagrams for a large sample 
of OH/IR stars Lewis et al. (2004) have recently sug¬ 
gested that the vast majority of them could have detached 
shells, possibly linked to a brief period of intense mass loss 
for oxygen-rich AGB stars (Lewis 2002). However, no de¬ 
tached CO shell has been found for any oxygen-rich object 
yet. 

Evidence of detached dust shells have been found in 
the cases of U Ant (Izumiura et al. 1997), U Hya (Waters 
et al. 1994), Y CVn (Izumiura et al. 1996), and R Hya 
(Hashimoto et al. 1998). All, but the last are carbon stars. 
Except for U Ant, there has been no detections of detached 
CO shells for these stars. 

Olofsson et al. (1990) suggested that the detached 
shells were an effect of a mass-loss-rate modulation caused 
by a He-shell flash (i.e., a thermal pulse). Vassiliadis & 
Wood (1993) studied this in some more detail, and fol¬ 
lowed the total mass-loss history of an AGB star. More 
recently, Schroder et al. (1998), Schroder et al. (1999) 
and Wachter et al. (2002) combined stellar evolutionary 
models and mass-loss prescriptions to study the effects on 
the evolution of an AGB star. They argued that during 
the low-mass-loss-rate evolution only carbon stars reach a 
critical (Eddington-like) luminosity that drives an intense 
mass ejection during a He-shell flash, hence explaining the 
absence (or at least rare occurrence) of detached shells to¬ 
wards M-type AGB stars. At higher mass-loss rates the 
detached shells become, in comparison, less conspicuous. 
Steffen & Schonberner (2000) used hydrodynamical simu¬ 
lations to follow the circumstellar evolution due to changes 
in the stellar mass-loss rate. They verified that a brief pe¬ 
riod of very high mass-loss rate translates into an expand¬ 
ing, geometrically thin shell around the star. They also 
concluded that an interacting-wind scenario is an equally 
viable explanation for the detached shells. 

Less dramatic modulations of the mass-loss rate on 
shorter time scales have been observed towards the 
archetypical high-mass-loss-rate carbon star IRG-l-10216 
(Mauron & Huggins 1999, 2000; Fong et al. 2003). Thin 
arcs, produced by mass-loss-rate modulations while the 
central star was ascending the AGB, are also seen around 
PPNe (Hrivnak et al. 2001; Su et al. 2003) and PNe 
(Gorradi et al. 2004). 

In this paper we present new single-dish millimetre line 
observations of R Scl, U Gam, V644 Sco, and DR Ser, 
Sect. 2. Of particular importance is the detection of a 
new carbon star with a young detached molecular shell, 
DR Ser, increasing the number of such stars to seven. In 
fact, there is not much hope of increasing this number 
further since most of the reasonably nearby AGB stars 
with mass loss have already been searched for circumstel¬ 
lar GO radio line emission. Section 3 describes the analy¬ 
sis of the molecular line and dust emission. Properties of 
the detached shells (Sect. 4), and the present-day winds 
(Sect. 5) of seven carbon stars with detached molecular 
shells. Table 1, are estimated. The results are discussed in 
Sect. 6, where different scenarios for producing the shells 
are reviewed. The usefulness of future high spatial reso- 


Table 1. Carbon stars with detached molecular shells. 


Source 

RA(2000.0) 
[h : m : s] 

DEC(2000.0) 

r : ' : "] 

Ref. 

R Scl 

01:26:58.09 

-32:32:35.5 

1 

U Cam 

03:41:48.17 

62:38:54.5 

2 

U Ant 

10:35:12.85 

-39:33:45.3 

3 

V644 Sco 

17:26:18.7 

-40:01:52 

4 

DR Ser 

18:47:21.02 

05:27:18.6 

3 

S Set 

18:50:20.04 

-07:54:27.4 

3 

TT Cyg 

19:40:57.00 

32:37:05.9 

5 


Refs. - 1) Wong et al. (2004), 2) Lindqvist et al. (1999), 3) 
Hipparcos catalogue, 4) Kholopov et al. (1985), 5) Olofsson 
et al. (2000). 


lution, as well as high-J GO observations, in further con¬ 
straining the shell characteristics are also discussed. Our 
conclusions are presented in Sect. 7. 

2. Observations 

2.1. Molecular line observations 

Olofsson et al. (1996) published molecular line spectra 
obtained towards all the stars in Table 1, but U Cam 
and DR Ser. Additional (sub)millimetre single-dish spec¬ 
tra used in the present analysis have been published in 
Olofsson et al. (1993a,b), Schoier & Olofsson (2001), and 
Wong et al. (2004). 

In this paper we present some complementary observa¬ 
tional data. The James Clerk Maxwell Telescope^ (JCMT) 
located at Mauna Kea was used to collect CO J = 3 ^ 2 
(345.796 GHz) and <7 = 4^ 3 (461.041 GHz) line emission 
from R Scl, DR Ser, and U Cam. The observations were 
performed in June, July, and October 2003, except for the 
J = 4 3 observations towards U Cam which were done 

in November 2000. CO J = 1 ^ 0 (115.271 GHz) and 
J = 2 ^ 1 (230.538 GHz) line emission from U Cam was 
obtained in November 1999 using the IRAM 30 m tele¬ 
scope^ at Pico Veleta. In addition, CO J = 1 ^ 0 line 
emission towards TT Cyg and ^^CO J = 2 ^ 1 (220.399 
GHz) line emission towards V644 Sco were detected using 
SEST in April 2002. The HCN J = 1 ^ 0 (88.632 GHz) 
data for DR Ser were collected during 1996-1997 using the 

^ The JCMT is operated by the Joint Astronomy Centre 
in Hilo, Hawaii on behalf of the present organizations: the 
Particle Physics and Astronomy Research Council in the 
United Kingdom, the National Research Council of Canada 
and the Netherlands Organization for Scientific Research. 

^ The IRAM 30 m telescope is operated by the Intitut de 
Radio Astronomie Millimtrique, which is supported by the 
Centre National de Recherche Scientiflque (France), the Max 
Planck Gesellschaft (Germany) and the Institute Geografico 
National (Spanin). 
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Onsala Space Observatory (OSO) 20m telescope^. HCN 
J = 1 ^ 0 and J = 3 ^ 2 (265.886 GHz) observations of 
V644 Sco were obtained using the Swedish-ESO submil¬ 
limetre telescope^ (SEST) at La Silla in April 2002. 

All observations were made in a dual beamswitch 
mode, where the source is alternately placed in the sig¬ 
nal and the reference beam, using a beam throw of about 
2' at JCMT, 4' at IRAM, 11' at OSO, and 12' at SEST. 
This method produces very flat baselines. Regular point¬ 
ing checks made on strong continuum sources (JCMT and 
IRAM) or SiO masers (OSO and SEST) were found to be 
consistent within wS". 

The observed spectra are presented in Fig. 1 and 
velocity-integrated intensities are reported in Table 2. The 
intensity scales are given in main beam brightness temper¬ 
ature, Tmb = T 4 /? 7 mbj where is the antenna tempera¬ 
ture corrected for atmospheric attenuation using the chop¬ 
per wheel method, and rymb is the main beam efficiency. 
For the CO J = 3 ^ 2 and J = 4 ^ 3 observations at 
the JCMT we have used main-beam efficiencies of 0.62 
and 0.5, respectively. For the IRAM CO J = 1 ^ 0 and 
J = 2 ^ 1 observations we have used main-beam effi¬ 
ciencies (i?eff/Feff) of 0.7 and 0.5, respectively. For the 
HCN J = 1 ^ 0 data obtained at OSO = 0-65 was 
used. For the SEST r/mb = 0.75 was used for the HCN 
J = 1 ^ 0 data, rymb = 0.7 for the CO J = 1 ^ 0 
data, and rymb = 0.5 for the HCN J = 3 ^ 2 and ^^CO 
J = 2 ^ 1 observations. The uncertainty in the absolute 
intensity scale is estimated to be about ±20%. 

The data was reduced in a standard way, by removing 
a first order baseline and then binned in order to improve 
the signal-to-noise ratio, using XS^. 

A summary of the available data, relevant to the anal¬ 
ysis in this paper, are given in Table 2. 

2.2. Continuum observations 

The observed spectral energy distribution of each source 
is used to constrain the dust modelling. We have selected 
observations that contain the total source flux at a spe¬ 
cific wavelength. In addition to IRAS fluxes (12-100 /im), 
near-infrared (NIR) JiLATM-photometric data have been 
obtained from the literature (see Kerschbaum 1999 and 
references therein). The IRAS and NIR data for a par¬ 
ticular source were generally not obtained at the same 
epoch, while the NIR photometry data was. No reliable 
100/im IRAS flux estimates are available for DR Ser and 
V644 Sco, due to cirrus contamination. For V644 Sco, 
DR Ser, S Set, and TT Cyg there exist observations from 


Table 2. Summary of single-dish millimetre line observations 
used in the analysis. 


Source 

Telescope 

Transition 

/ Tmbdv 

[Kkms“^] 

Ref. 

R Scl 

IRAM 

CO(l -0) 

76.2 

4 


IRAM 

CO(2 - 1) 

67.7 

4 


SEST 

CO(l-O) 

26.1 

4 


SEST 

CO(2- 1) 

47.8 

4 


SEST 

CO(3-2) 

63.7 

3 


JCMT 

CO(2- 1) 

50.6 

6 


JCMT 

CO(3 - 2) 

62.8 

1 


JCMT 

CO(4-3) 

63.2 

1 


SEST 

i®CO(l -0) 

1.9 

5 


SEST 

i®CO(2- 1) 

4.2 

5 


SEST 

HCN(1 - 0) 

1.3 

3 


SEST 

HCN(3 - 2) 

7.7 

3 


SEST 

HCN(4 - 3) 

3.9 

3 


HHT 

HCN(3- 2) 

2.4 

7 


HHT 

HCN (4 - 3) 

4.1 

7 

U Cam 

IRAM 

CO(l-O) 

21.8 

1 


IRAM 

CO(2- 1) 

97.9 

1 


JCMT 

CO(3-2) 

66.7 

1 


JCMT 

CO(4-3) 

54.6 

1 

U Ant 

SEST 

CO(l-O) 

10.8 

3 


SEST 

CO(2- 1) 

13.1 

3 


SEST 

CO(3-2) 

9.7 

3 

V644 Sco 

SEST 

CO(l-O) 

7.6 

3 


SEST 

CO(2- 1) 

20.7 

3 


SEST 

CO(3-2) 

21.4 

3 


JCMT 

CO(2- 1) 

22.2 

6 


JCMT 

i®CO(2- 1) 

1.1 

1 


SEST 

HCN(1 - 0) 

0.3 

1 


SEST 

HCN(3 - 2) 

0.4 

1 

DR Ser 

SEST 

CO(l-O) 

3.7 

2 


NRAO 

CO(2- 1) 

11.3 

2 


JCMT 

CO(3-2) 

14.7 

1 


JCMT 

CO(4-3) 

8.5 

1 


OSO 

HCN(1 - 0) 

0.6 

1 

S Set 

SEST 

CO(l-O) 

5.8 

3 


SEST 

CO(2- 1) 

4.3 

3 


SEST 

CO(3 - 2) 

3.1 

3 


SEST 

i3CO(l -0) 

0.2 

3 

TT Cyg 

SEST 

CO(l-O) 

5.0 

1 


OSO 

CO(l-O) 

6.4 

3 


IRAM 

CO(l-O) 

7.5 

3 


IRAM 

CO(2- 1) 

6.0 

3 


Refs. - 1) This paper. 2) Schoier & Olofsson (2001). 3) Olofsson 
et al. (1996). 4) Olofsson et al. (1993b). 5) Schoier & Olofsson 
(2000). 6) JCMT public archive. 7) Sieging (2001). 


® The Onsala 20 m telescope is operated by the Swedish 

National Facility for Radio Astronomy, Onsala Space obser- , i i o ^ /iv 

. . r , \ the Midcourse Space Experiment (MSX) in the range 

vatory at Chalmers University of technology or.i i ■ 

4 rni oT^orn i j • • o j• 1 AT i• 1 o ~ 11 uTo. wFich Lave been used in the analysis. 

The SEST was operated jointly by the Swedish National ^ 

Facility for Radio Astronomy and the European Southern ^ Set has the largest detached shell of the sample 
Observatory (ESO). sources with a diameter of approximately 140". This is 

® XS is a package developed by P. Bergman to reduce and comparable to the size of the IRAS beam at 60 /im 
analyse a large number of single-dish spectra. It is publically (120") and resolution effects start to become important, 
available from ftp://yggdrasil.oso.ch.almers.se Groenewegen & de Jong (1994), in their modelling of 
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Fig. 1. New observations of CO and HCN radio line emission. 
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Fig. 2. Observed SEDs for the seven carbon stars with detected detached molecular shells. Best-fit models to the observed 
SEDs are overlayed (solid line). The model SEDs are composed of a stellar part (dashed line), a present-day wind (dotted line), 
and a detached shell (dash-dotted line). The model SEDs have been corrected for interstellar extinction. Stellar and envelope 
parameters are summarized in Tables 3 and 4. 


S Set, estimated that about 30% of the flux was resolved 
out. Given the larger beam (220") the effect at 100 /im 
is much lower as it is for all other sources in the sample. 
For this reason resolution effects has not been taken into 
account in the analysis. 

We have obtained archival data taken by the Short 
Wavelength Spectrometer (SWS) onboard the Infrared 
Space Observatory (ISO) for U Cam and estimated the 
flux at 40 ^m. U Cam is the source with the smallest 
detached shell of 15" in diameter, smaller than the ISO 
aperture (17" x 40") at 40 ^m. 

The observed SEDs are presented in Fig. 2. 


3. Radiative transfer 


ters describing the shell need to be varied together with 
the mass-loss rate describing the present-day wind. 

In what follows the procedures applied to determine 
the properties of the shells, as well as those of the present- 
day winds, from both continuum and millimetre line ob¬ 
servations, are descibed. 

The best fit model is found by minimizing 


N 

i=l 

where F\ is the flux, and a\ the uncertainty in the mea¬ 
sured flux, at wavelength A, and the summation is done 
over N independent observations. 


(Tmod.A Eobs,A) 


* Z 




( 1 ) 


This section describes the basic assumptions made in the 
analysis of the observed molecular line and continuum 
emission, and the methods used in the treatment of the 
radiative transfer. 

3.1. General considerations 

The CSEs are assumed to consist of two spherically sym¬ 
metric components: an attached CSE (aCSE), i.e., a nor¬ 
mal continuous AGB wind, and a detached, geometrically 
thin, CSE (dCSE) of constant density and temperature. 
Both components are assumed to be expanding at con¬ 
stant, but different, velocities. From the observed line 
spectra and SEDs it is clear that in many cases it is not 
easy to separate the contribution from the present-day 
wind and the detached shell. In these cases the parame- 


3.2. Molecular line modelling 

A detailed non-LTE radiative transfer code, based on the 
Monte Carlo method (Bernes 1979), is used to perform the 
excitation analysis and to model the observed circumstel- 
lar line emission. The code is described in detail in Schoier 
& Olofsson (2001) and has been benchmarked, to high ac¬ 
curacy, against a wide variety of molecular line radiative 
transfer codes in van Zadelhoff et al. (2002) and van der 
Tak et al. (in prep.). 

The excitation analysis includes radiative excitation 
through vibrationally excited states, and a full treat¬ 
ment of line overlaps between various hyperflne com¬ 
ponents in the case of HCN, see e.g. Lindqvist et al. 
(2000). Relevant molecular data are summarized in 
Schoier et al. (2005) and are made publicly available at 
http: //www. strw. leidenuniv .nl/^moldata. 
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Table 3. Stellar properties. 


Source 

Var. type 

P 

[days] 

D 

[pc] 

L* 

[L©1 

T* 

[K] 

R Scl 

SRb 

370 

290^ 

4300 

2625 

U Cam 

SRb 

400 

430^ 

7000 

2695 

U Ant 

Lb 

— 

260^ 

5800 

2810 

V644 Sco 

Lb 

— 

700® 

4000 

2615 

DR Ser 

Lb 

— 

760® 

4000 

2650 

S Set 

SRb 

148 

400® 

4900 

2755 

TT Cyg 

SRb 

118 

610® 

3200 

2825 


^ Hipparcos distance (Knapp et al. 2003) 

^ From period-luminosity relation (Knapp et al. 2003) 
® From adopting L* = 4000 L© 


3.3. Continuum modelling 

An independent estimate of the properties of the detached 
shell and the present-day wind can be made from the ex¬ 
cess emission observed at wavelengths longer than « 5 /rm 
and 20 gm, respectively. The dust emissions from both 
the shell and the present-day wind are expected to be op¬ 
tically thin at all relevant wavelengths (> 1 /xm). The as¬ 
sumption of optically thin emission greatly simplifies the 
radiative transfer of the continuum emission. The simple 
model presented in detail below has been verihed against 
a full radiative transfer code (Schoier, in prep.) for the 
densest winds encountered in this project, and they are 
found to agree within the numerical accuracy of the code 
itself (< 5%). 

In the optically thin limit the contribution to the 
SED can be separated into a stellar component a 

component from the present-day wind and that 

from the detached shell which are solved for in¬ 

dividually. The total SED is then obtained from F^, = 
^star _|_ ^wind _|_ ^sheii violation of energy 

conservation introduced in this way is negligible, ^5%). 
Note that we ignore any beam hlling effects (see Sect. 2.2). 

The flux received by an observer located at a dis¬ 

tance D from an expanding, spherically symmetric wind 
produced by a constant rate of mass-loss is 

- J ndust{r)B„{Tdust{r))r^dr, (2) 

where Qi, is the absorption efficiency of a spherical dust 
grain with radius Ugr, n-dust the number density of dust 
grains, the blackbody brightness at the dust temper¬ 
ature Tdust, and the integration is made from the inner 
(Ri) to the outer (Re) radius of the CSE. 

Eor an extended, geometrically thin, isothermal shell 
at constant density, adequate for describing the detached 
shells treated here, Eq. (2) reduces to 

j^shell ^ ^^'^^, F.(rdust)Mdust, (3) 


where the dust mass of the shell, Mdust, and the specific 
density of a grain, pgr, have been introduced. 

The local dust temperature, Tdust (?'), is determined 
by balancing the heating (iJ) provided by absorption 
of infalling radiation and cooling (C) by subsequent re¬ 
emission. The heating term can be expressed as 

H = TTal^ridustir) J Fj;^^^{r)Q^du, (4) 

where the dust grains are assumed to be of the same size 
and spherical. The dust grains are assumed to locally radi¬ 
ate thermal emission according to Kirchoff’s law and the 
resulting cooling term is 

C = 47r^agrndust(?') J B^{Tdust ir))Q„dv. (5) 


Modelling the SED provides information only on the 
density structure of the dust grains, oc Mdust/^’dust, where 
Md is the dust-mass-loss rate and Wdust the dust veloc¬ 
ity. In order to calculate Mdust, and eventually the gas 
mass-loss rate Mgas, one has to rely upon some model de¬ 
scribing the dynamics of the wind. There are now growing 
evidence that the winds of AGB stars are driven by radia¬ 
tion pressure from stellar photons exerted on dust grains. 
Through momentum-coupling between the dust and gas, 
molecules are dragged along by the grains. The stationary 
solution to such a problem has been extensively studied 
(e.g., Habing et al. 1994). The terminal (at large radii) gas 
velocity, Vgas, is given by 

vgus = V2GM*(r - l)/ro, (6) 


where G is the gravitational constant, M* the stellar mass, 
T* the luminosity, and T the ratio of the drag force and 
the gravitational force on the gas, given by 

Y _ 3(QF/agr)T^,4t Ugas .js 

IbTrGcMpgr Udust’ 

where Qf is the flux-averaged momentum-transfer efh- 
ciency, c the speed of light, and the dust-to-gas mass- 
loss-rate ratio. The drift velocity x;dr = ^’dust — Vgas can be 
obtained from 


Vdr = 




CM„; 


( 8 ) 


In the modelling, amorphous carbon dust grains 
with the optical constants presented in Suh (2000) were 
adopted. The properties of these grains were shown to 
reproduce, reasonably well, the SEDs from a sample of 
carbon stars. Eor simplicity, the dust grains are assumed 
to be spherical and of the same size with Ogr = 0.1 pm and 
Pgr = 2gcm“^. The absorption and scattering efficiencies 
were then calculated using standard Mie theory (Bohren 
& Huffman 1983). Eor the dynamical solution we further 
assume a stellar mass of M = 1 Mq and adopt an inner 
radius of the envelope tq = 8x 10^^ cm (corresponding to a 
dust condensation temperature of « 1200K). 
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4. Dust modelling 

The state of the dust grains will have a great impact on the 
gas modelling, affecting both its dynamics and molecular 
excitation. However, the reverse is not true so the dust 
modelling can safely be done without bothering about the 
state of the gas, except for the terminal gas velocity which 
enters in the solution of the dynamics and allows the total 
mass-loss rate to be estimated. Fortunately, the terminal 
gas velocity is an observable. 

4.1. Stellar parameters 

The first step in the modelling is to determine the lumi¬ 
nosity (alternatively, the distance) from htting the SED 
shortward of w 10 pm where the stellar emission fully dom¬ 
inates in these optically thin envelopes. The stellar emis¬ 
sion, is assumed to be that of a blackbody. The stel¬ 
lar effective temperatures, T*, for our sample sources were 
adopted from Bergeat et al. (2001). Accurate Hipparcos 
parallaxes exist only for U Ant, S Set, and TT Cyg (Knapp 
et al. 2003), and other means of determining the distance 
to the remaining objects need to be adopted. For R Scl and 
U Cam we used the period-luminosity relation by Knapp 
et al. (2003) to obtain the distances, adopting the peri¬ 
ods listed by Kholopov et al. (1985) and Knapp et al. 
(2003). The distances to the irregular variables V644 Sco 
and DR Ser were obtained after adopting a luminosity of 
4000 Lq. All the stellar parameters are listed in Table 3. In 
the htting procedure, interstellar extinction was accounted 
for by adopting the extinction law of Cardelli et al. (1989) 
and using the vales for extinction in the visual. Ay, listed 
in Bergeat et al. (2001). 

4.2. Detached shell and present-day wind 

The observational constraints, in the form of SEDs cov¬ 
ering the wavelength range w 10 — 100 pm are analyzed 
using the y^-statistic dehned in Eq. (1). A calibration 
uncertainty of 20% was added to the total error budget 
and will, in most cases, dominate the error which in turn 
means that huxes at various wavelengths typically have 
the same weight in the y^-analysis. It should be noted 
that the present-day mass loss is more difficult to esti¬ 
mate than the properties of the detached shell given its 
lower contrast with respect to the stellar contribution. 

For the present-day wind there is only one adjustable 
parameter in the modelling, the gas-mass-loss rate. The 
emission from the shell is fully described by adjusting the 
temperature of the dust grains (Tdust) and the total dust 
mass in the shell (Mdust)- The best fit models all have re¬ 
duced values of « 1, indicating good hts. The param¬ 
eters obtained are listed in Table 4. The dust-to-gas ratio 
(4^) in the present-day wind obtained from the dynamical 
model is also listed. The sensitivity of the results to the ad¬ 
justable parameters are illustrated in Fig. 3 for the young 
shell around DR Ser and the older shell around TT Cyg. 
The final hts to the observed SEDs are presented in Fig. 2, 


where also the contribution of each of the three compo¬ 
nents (stellar, present-day wind, and detached shell) are 
indicated. 

The dust temperature in the shell can be translated 
into a radial distance of the shell from balancing the 
heating [Eq. (4)], provided by the central star, and the 
cooling [Eq. (5)], provided by the re-emission of the ob¬ 
served photons at longer wavelengths, as in the case of 
the present-day wind. The distances obtained are given in 
Table 4. As shown in Fig. 4, although the errors involved 
are quite large, there is a clear correlation between the ra¬ 
dial distance of the shell and its mass, Mdust fx 
Pearson’s correlation coefficient r is 0.97 indicating an al¬ 
most complete positive correlation. The fact that the mass 
increases with distance suggests that matter is being swept 
up as the shell expands away from the star. Such a sce¬ 
nario will be further elaborated in Sect. 6.1. 

U Cam is the source which has the lowest contrast 
between the stellar and shell contributions to the SED. 
It is the youngest shell in the sample located close to the 
central star. Sect. 5., complicating also its separation from 
the present day wind. This makes it hard from the dust 
modelling to constrain the properties of the circumstellar 
dust for this particular source. 

Groenewegen & de Jong (1994) modelled the SED of 
S Set and found a shell dust mass of « 2x lO”"* M©. Given 
the uncertainties in the adopted dust properties, this is 
consistent with the mass derived in the present analysis. 
Placing the dust shell at the same location as that of CO 
we derive a dust mass of «8xl0“® Mq for our somewhat 
smaller distance. 

Izumiura et al. (1997) in their analysis of HIRAS im¬ 
ages found evidence for two separate dust shells around 
U Ant where the inner of the two shells is coinciding with 
that of CO. They derive a total mass for this shell of 
« 5x 10“^ M© consistent within the uncertainties with our 
estimate of «2xl0“^ M© from the CO modelling. The 
outer shell, which has no known molecular counterpart, is 
similarly thought to have been produced in another ther¬ 
mal pulse some ~ 10^ yr before. The dust mass contained 
in this shell is estimated to be «4xl0“^ M© and its con¬ 
tribution to the SED is «15% at 60 pm. 

5. Molecular line modelling 

5.1. Line profiles 

Many of the observed CO spectra exhibit a triple-peaked 
line profile, e.g., Olofsson et al. (1996). The, usually nar¬ 
row, peaks at the extreme velocities arise from a spatially 
resolved shell, whereas the central component is produced 
by the present-day wind, i.e., a normal low-mass-loss-rate 
AGB wind. Since the present-day wind is in general ex¬ 
panding at a significantly lower velocity, Wp, than that of 
the shell, 14, emission at velocities between u = u* ± Up 
and V = V* ± 14 should contain only the contribution from 
the shell. Each of these velocity intervals is further divided 
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Table 4. Results from dust modelling (with Ict errors) 


dCSE aCSE 


Source Mdust Tdust Ra Mgaa t!' 

[Mq] [K] [cm] [Mgyr-I] 


R Scl 

(3.2 ± 

2.0) 

X 

10“® 

75 

± 

15 

(1.2 ± 

0.6) 

X 

lO^" 


3.8 

X 

10“ 

-7 

1.7 

X 10"® 

U Cam 

1.5 X 10“ 

6 _ 

l.C 

) X 10"® 

70 

± 

45 


— 




4.0 

X 

10“ 

-7 

1.2 

X 10“® 

U Ant 

(1.3 ± 

1.2) 

X 

lO"-* 

55 

± 

20 

(5.1 ± 

4.0) 

X 

lO^" 


7.0 

X 

10“ 

-8 

3.0 

X 10~® 

V644 Sco 

(1.4 ± 

0.9) 

X 

lO""* 

61 

± 

9 

(1.8 ± 

0.7) 

X 

10^'^ 

< 

5.0 

X 

10' 

-8 


— 

DR Ser 

(3.5 ± 

2.5) 

X 

10"® 

68 

± 

13 

(1.5 ± 

0.7) 

X 

lOi" 

< 

7.0 

X 

10' 

-8 


— 

S Set 

(8.0 ± 

7.0) 

X 

lO""* 

34 

± 

9 

(1.3 ± 

0.8) 

X 

lO^® 

< 

8.0 

X 

10' 

-8 


— 

TT Cyg 

(4.3 ± 

3.7) 

X 

lO"-* 

39 

± 

10 

(6.9 ± 

4.1) 

X 

lOi" 

< 

3.2 

X 

10' 

-8 


— 




Fig. 3. -maps showing the sensitivity of the adjustable pa¬ 
rameters Tdust and Mdust in the dust modelling of the detached 
shells around DR Ser and TT Cyg. Contours are drawn at 
Xmin + (2.3,6.2,11.8) indicating the 68% (Ter’), 95% (‘2cr’), 
and 99.7% (‘3cr’) confidence levels, respectively. The quality 
of the best fit model can be estimated from the reduced chi- 
squared statistic Xred=Xmin/(.^~2) shown in the upper right 
corner. 


into two separate components so that the total x^"Statistic 
also contains information on the line prohle. 



10 ^^ 10 ^® 
Rg [cm] 


Fig. 4. The shell dust (Mdust; upper panel) and gas (Mgas; 
lower panel) masses plotted as a function of the radius of the 
shell {Ra). The solid line shows the mass-radius relation of the 
interacting wind model presented in Sect. 6.2. 


Information on the properties of the present-day mass- 
loss rate is contained in the intensity in the range v = 
ti* ± rip. However, there is also generally a contribution 
here from the detached shell and a good model for the 
shell is required in order to investigate the present-day 
wind. This makes it particularly difficult to separate the 
two components in all cases except U Ant, S Set, and 
TT Cyg, (e.g., Olofsson et al. 1996). In the case of U Cam 
the present-day wind is separated in the CO line interfer- 
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ometer data (Lindqvist et al. 1999). In the cases of R Scl, 
V644 Sco, and DR Ser the detected HCN emissions indi¬ 
cate expansion velocities about a factor of two lower than 
in the case of the CO emission. The HCN emission is inter¬ 
preted as arising predominantly in the present-day wind. 
This was also the conclusion reached by Wong et al. (2004) 
for R Scl. 

5.2. The detached shells 

In the analysis there exist four adjustable parameters: the 
H 2 number density (riHa), kinetic gas temperature (Tkin), 
thickness (ARs) and location of the shell (Rs). The frac¬ 
tional abundance of CO (in relation to H 2 ), /co, is as¬ 
sumed to be equal to 1.0x10“^. In three cases (TT Cyg, 
Olofsson et al. (2000); U Cam, Lindqvist et al. (1999); S 
Set, Olofsson et ah, in prep.) the sizes of the shells are well 
known from interferometric CO observations. In the case 
of U Ant the size of the shell is reasonably well known from 
single-dish CO mapping. R Scl has a prominent detached 
shell at « 8.7x10^® cm seen in scattered light (Gonzalez 
Delgado et al. 2001). In these five cases the position of 
the shell was fixed to match the observations. Although 
high-spatial-resolution observations exist, little is known 
about the actual thickness of these shells. CO interfer¬ 
ometric observations do not resolve the radial thickness 
of the shells. Here we adopt a value of 1.0x10^® cm. The 
best-fit model is found from minimizing the y^-statistic as 
described above. 

As an example, in Fig. 5 y^-maps obtained from the 
modelling of the shell around V644 Sco are shown. The 
location of the shell (Rs) is unknown and used as an ad¬ 
justable parameter. However, it turns out that the shell ra¬ 
dial distance from the central star is the best constrained 
parameter in the modelling, together with the mass con¬ 
tained in the shell. Fig. 5 (lower panel) shows the sensitiv¬ 
ity of the model to changes in Rs and nua for a fixed tem¬ 
perature of 100 K. In a broad temperature range we find 
that good fits are obtained for i?s=(l-l±0.2)xl0^^cm. 
Fixing the radial distance to 1.1x10^^ cm, the density 
and temperature are constrained to 1050 ± 150 cm“® and 
170± 100 K, respectively (Fig. 5; upper panel). The mass 
contained in the shell is Mgas = (2.5±0.3) x 10“® Mq. 

The derived mass in the shell is not very sensitive to 
the adopted thickness of the shell since the emission in 
almost all observed lines is optically thin. The density in 
the shell, however, goes roughly as oc There is 

also a weak temperature dependence. 

The CO line emission from the lower rotational transi¬ 
tions is not very sensitive to changes in the kinetic gas 
temperature once Tkin > 50 — lOOK. In order to better 
constrain the temperature in the shell transitions involv¬ 
ing higher J-levels needs to be observed as illustrated in 
Sect. 6.4. 

The properties of all dCSEs determined from the CO 
line modelling are summarized in Table 5. Fig. 4 illustrates 
that, as in the case of the dust emission, there is a clear 


Ms [Mol 

10“® 2x10“^ 5x10“^ 




Fig. 5. y^-maps showing the sensitivity of the adjustable pa¬ 
rameters riHa, Tkin, and As in the CO modelling of the de¬ 
tached shell around V644 Sco. Contours are drawn at y^in + 
(2.3, 6.2, 11.8) indicating the 68% (Ter’), 95% (‘2a’), and 99.7% 
(‘Sc’) confidence levels, respectively. The quality of the best fit 
model can be estimated from the reduced chi-squared statistic 
Xred=Xmin/(N—2) shown in the lower left corner. Also indi¬ 
cated is the value of the parameters that is kept hxed in these 
maps. 

trend that the gas mass of the shell increases with its radial 
distance from the central star. Again a strong correlation 
is found Mgas oc i^O-SSio.is _ 0.93) further supporting 
a scenario where matter is being swept up as discussed 
in Sect. 6.1. Also in indicated in Table 5 is the age of 
the dCSEs estimated from Rs/Vg. Since there is a clear 
trend that 14 decrease with Rg (Sect. 6.1) these ages are 
strictly upper limits to the age since the expansion was 
likely faster at earlier epochs. 

5.3. The present-day winds 

Once the properties of the detached shell are determined 
and its contribution to the observed CO line emission is 
estimated, the corresponding estimates for the present- 
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Table 5. Results from molecular line modelling (with Ict errors bars) 


Source 


dCSE 




aCSE 


Mgas 

[M©] 

2kin 

[K] 

Rs 

[cm] 

Fs 

[kms"®] 

age 

[yr] 

Mgas 

[M0yr"®] 

Vp 

[kms"®] 

Molecule® 

R Scl 

(2.5 ±0.5) X 10'® 

> 20 

8.7 X 10®® 

15.5 

1800 

3.0 X 10"’® 

10.5 

HCN 

U Cam 

(9.5 ± 1.0) X 10-^ 

> 130 

4.7 X 10®® 

23.0 

650 

2.0 X 10"'® 

12.0 

CO 

U Ant 

(1.9 ±0.2) X 10"® 

> 200 

1.7 X 10®^ 

19.0 

2800 

2.0 X 10"® 

4.0 

CO 

V644 Sco 

(2.5 ±0.3) X 10"® 

170 ± 100 

(1.1 ±0.2) X lO®®" 

23.0 

1500 

5.0 X 10"® 

5.0 

HCN 

DR Ser 

(1.0 ±0.2) X 10"® 

> 100 

(8.0 ± 1.5) X 10®® 

20.0 

1300 

3.0 X 10"® 

5.0 

HCN 

S Set 

(7.5 ± 1.0) X 10"® 

> 60 

4.2 X 10®®" 

16.5 

8100 

2.0 X 10"® 

4.0 

CO 

TT Cyg 

(4.0 ± 1.0) X 10"® 

> 200 

3.2 X 10®®" 

12.5 

8100 

3.2 X 10"® 

4.0 

CO 


^ The molecular line emission used for estimating the properties of the aCSEs. Note that in the analysis of the dCSEs CO line 
emission was used in all cases. 





-50 0 50 






Fig. 6. Best fit models (solid line) including both the present-day wind and detached shell emission overlayed on observed 
spectra (histograms) for U Ant and U Cam. The uncertainty in the intensity scale is ±20%. The narrow spikes found near the 
blue wing in the spectra of U Cam are of interstellar origin. 


day wind can be can be made. In the modelling of the 
present-day wind the temperature structure is solved for 
self-consistently as described in Schoier & Olofsson (2001) 
using the results from the dust modelling performed in 
Sect. 4 to describe the properties of the dust entering in 
the heating of the gas envelope. It should be emphasised 
that the actual CO line intensities are not very sensitive 
to the temperature structure since the excitation is domi¬ 
nated by the radiation field from the central star in these 
low-mass-loss-rate objects (Schoier & Olofsson 2001). 

Adopting a CO abundance of 1.0x10“^ relative to H 2 
leaves just the mass-loss rate, Mgas, as an adjustable pa¬ 
rameter in the modelling. The expansion velocity of the 
present-day wind is fixed from matching the observed lines 
with the model (assuming a micro-turbulent line broad¬ 
ening of 1.0 kms“^, in addition to the thermal motion). 
The CO envelope size is calculated from the photodissoci¬ 


ation model of Mamon et al. (1988). The applicability of 
these results to the present-day wind in these objects is a 
bit questionable since the shell itself may provide shield¬ 
ing against photodissociation. In the modelling, the CO 
J = 1 ^ 0 intensity is most sensitive to the size of the 
envelope (Schoier & Olofsson 2001). 

For some of the younger detached shell sources like 
R Scl, V644 Sco, and DR Ser it is difficult, given the cur¬ 
rent set of CO observations available, to estimate reliable 
mass-loss rates. In these objects the observed HCN line 
emission has been used instead. HCN is only probing the 
present-day wind since it is readily photodissociated at the 
distances where the detached shells are found. Note, how¬ 
ever, that there is some indication that HCN J = 1 —> 0 
emission from the shell is detected in V644 Sco (Fig. 1). 
As in the case of the CO modelling an assumption of the 
HCN abundance distribution needs to be made in order to 
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estimate the mass-loss rate. HCN is a species with a pho- 
tospheric origin. Olofsson et al. (1993a) estimated HCN 
abundances in the photospheres of four of the sample stars 
and obtained values of ~ 0.8 — 5.4 x 10“® relative to H 2 . 
Here we adopt a generic value of 2 x 10“® for the anal¬ 
ysis. The size of the HCN envelope was then calculated 
from a photodissociation model based on the results from 
Lindqvist et al. (2000). It should be pointed out that for 
each HCN model a CO model using the same mass-loss 
rate needs to be calculated first in order to obtain the cor¬ 
rect temperature profile to use. This also makes it possible 
to check if the CO model is consistent with the observed 
CO lines. In all three cases it is. 

U Cam is the only source where the abundance of HCN 
in the aCSE can be determined based on the results from 
the CO modelling. Using the mass-loss rate and expan¬ 
sion velocities from Table 5, an HCN abundance of ap¬ 
proximately 4 X 10“® (relative to that of H 2 ) is obtained 
using the results from Lindqvist et al. (2000) to calculate 
the HCN envelope size. A more detailed investigation of 
the HCN and CN envelopes around U Cam based on new 
high-spatial-resolution interferometric observations is un¬ 
derway (Lindqvist et al., in prep.). U Cam, together with 
V644 Sco, are the only sources where HCN emission from 
the detached shell is detected. 

The expansion velocities and derived mass-loss rates 
for the present-day wind (aCSE) are presented in Table 5. 
The uncertainties in the mass-loss-rate estimates, within 
the adopted circumstellar model, are about ±50% for es¬ 
timates based on CO observations and as high as a factor 
~ 5 when HCN emission is used, mainly due to the larger 
uncertainty in the adopted photospheric HCN abundance. 
Examples of line profiles taken from the best-fit models, 
including both contribution from the present day wind and 
detached shell, are given in Fig. 6 for U Cam and U Ant. 

In Fig. 7 the present-day-wind mass-loss rates and 
gas expansion velocities are plotted against the size of 
the detached shell. Even though the uncertainties in the 
mass-loss rates are substantial (particularly those based 
on HCN emission) we conclude that there is a trend of 
decreasing present-day mass-loss rate with increasing size 
of the shell. This is further corroborated by the similar 
decrease in the gas expansion velocity with increasing size 
of the shell. The velocity is an observational property that 
can be determined with high accuracy, and it is known to 
correlate reasonably well with the mass-loss rate for ‘nor¬ 
mal’ CSEs (e.g., Olofsson 2003). The implications of this 
result is further discussed in Sect. 6. 

5.4. The detached shell around R Scl 

R Scl deserves a separate discussion. From the dust mod¬ 
elling it appears to be a normal detached shell source 
(see Table 4). However, the CO line modelling of the 
main isotope places the shell at about half the distance 
of what is actually determined from scattered-light ob¬ 
servations (Gonzalez Delgado et al. 2001). The mass of 



2x10^® 5x10^® 10^^ 2x10^'^ 5x10^'^ 

Rg [cm] 


Fig. 7. The present-day-wind gas expansion velocity (upper 
panel) and mass-loss rate (lower panel) plotted as a function 
of the radius of the shell (Rs). 

the shell is also unusually high « 1.0x10 “^Mq for such 
a young shell. Modelling of the rarer ^®CO isotopomer, 
however, gives results consistent with the larger radial 
position of the shell, and a much lower shell mass of ap¬ 
proximately 2.5xlO“^M0 is obtained [adopting the pho¬ 
tospheric ^^C/^^C ratio of 19 determined by Lambert et al. 
(1986)]. The properties of the detached shell around R Scl 
reported in Table 5 are based on the ^^CO modelling. 

The reason for the discrepancy obtained from the ^^CO 
modelling is not known. Possibly, the mass-loss-rate his¬ 
tory of this object is more complicated. Therefore, detailed 
interferometric CO observations are important for resolv¬ 
ing this issue. 

6. Discussion 

6.1. Evidence for interacting winds 

As shown in Fig. 4 there is compelling evidence from both 
the dust and gas modelling that a detached shell gains 
mass as it expands away from the central star. This would 
indicate that material is being swept up in the process. 
Further support for such a scenario comes from the ob- 
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Fig. 8. The shell expansion velocity (14) plotted against the 
shell radius {Rb', as determined from CO line data). The solid 
line shows the expected velocity-radius relation of the interact¬ 
ing wind model presented in Sect. 6.2 



Rs [cm] 

Fig. 9. The dust-to-gas ratio (T) plotted against the radius of 
the shell {Rs', assuming that the dust and gas shells coincide 
spatially). 

served width of the CO line emission emanating from the 
shell. Fig. 8 shows a correlation (r = —0.64) between the 
expansion velocity of the shell, 14, and its location, Rs 
in that the farther away the shell is located the lower its 
velocity. This could be naturally explained if the shell is 
indeed interacting with some other medium. In addition, 
the high inferred kinetic temperatures (actually lower lim¬ 
its) are also a signpost of interaction. 

The strongest argument for the Ms-Rs relation is 
that it is obtained in two completely independent ways. 
Nevertheless, the possibilities of systematic errors in the 
analysis and observational biases should be investigated. 
One way of looking for possible systematic errors in the 
analysis, which also includes both the dust and the gas es¬ 
timates, is to calculate the dust-to-gas mass-loss-rate ratio 



Rs [cm] (gas) 

Fig. 10. Relation between the radius of the shell as determined 
from dust and gas modelling, respectively. 

(dt) in the shells. To do this we have placed both the gas 
and dust shells at the same radial distance. The resulting 
values for dt as a function of the radial distance of the 
shell are shown in Fig. 9. There is no apparent correla¬ 
tion (r = 0.49) between the two quantities, and the results 
are consistent with the same value of T in all the shells, 
ss 0.012. This result provides no indication of a systematic 
error. On the other hand, the dust-to-gas mass-loss-rate 
ratios derived for the present day winds (Table 4) are much 
lower (by a factor 3-10) than those found for the shells, 
but fully consistent with those normally quoted for low 
to intermediate mass-loss-rate carbon stars (Groenewegen 
et al. 1998). As shown by Groenewegen et al. (1998) there 
is a trend that the dust-to-gas ratio increases with mass- 
loss rate and that for M > 5x 10“® Mq yr“^ (representa¬ 
tive of the high-mass-loss-rate epoch producing the de¬ 
tached shell) T* ~ 0.01. 

The fact that the locations of the dGSEs as determined 
from the dust modelling correlate very well with those esti¬ 
mated from GO observations, as shown in Fig. 10, suggests 
that the dust and gas shells are located relatively close 
spatially. Gonzalez Delgado et al. (2003) in their study of 
scattered stellar light in the GSE around U Ant found that 
the dust shell is located about 7" further out than that of 
GO (corresponding to a drift velocity of « 3kms“^). 

The most important parameter in the dust analysis 
is the dust temperature. The derived dust mass scales as 
TJ]ist in fiis Rayleigh-Jeans limit and as exp{hi'/kT^ust) in 
the Wien’s regime, and the inferred shell radius as 
(if the grain emissivity scales as u). This means that even 
in the case of identical shells a spread in estimated dust 
temperatures will result in a mass-radius relation (with a 
scatter depending on the properties of the central stars). 
Thus, in the dust analysis the accuracy of the dust temper¬ 
ature estimate is crucial. As shown in Table 4 these are 
rather well constrained (even though longer-wavelength 
data, free of cirrus contamination, are highly desirable), 
and we conclude that the dust temperature estimates are 
not the reason for the mass-radius relation. 
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There is possibly an observational bias in the sense 
that low-mass shells of large radii would be difficult to de¬ 
tect in both dust and CO line emission. We cannot exclude 
that such shells exist. On the other hand, we see no reason 
(other than the scarcity of stars with detached shells) why 
high-mass shells of small size should be missed, and are 
therefore reasonably certain about the upper envelope of 
the shell mass for a given shell size. 

6.2. A simple interacting-wind scenario 

Thus, we conclude that the observed mass-radius relations 
are most likely correct, and use a simple interacting-wind 
model to estimate whether our data can be explained in 
such a scenario using reasonable assumptions on AGB 
winds. Kwok et al. (1978) presented the results of a fully 
momentum-coupled wind interaction. In the model a fast 
moving wind is running into a slower moving wind. Both 
winds are assumed to be spherically symmetric, have con¬ 
stant mass-loss rates, and expand at constant velocities. 
The qualitative validity of this simple picture has been 
verified by detailed modelling performed by Steffen & 
Schonberner (2000). 

In our case it should be noted that the high velocities 
observed for the younger detached shell sources strongly 
suggests that the mass-loss rate forming the fast moving 
wind was also very high, say about 10 “® M 0 yr“^, since 
the mass-loss rate and gas expansion velocity correlates 
positively for normal CSEs (e.g., Olofsson 2003). Based 
on Fig. 4 it is reasonable to assume that the initial mass 
of the shell is of the order 10“^ Mq, i.e., with a mass- 
loss rate of 10 “^ Mq yr“^ the ejection period is about 10 ^ 
years. Thus, a brief period of very high mass-loss rate is 
inferred producing a fast moving shell. 

We have modihed the model of Kwok et al. (1978) to 
treat a scenario in which a fast moving shell is running 
into a slower moving wind. Conservation of the total mo¬ 
mentum contained in the shell (both ejected and swept-up 
material) gives an expression for the shell velocity 14 

^ MejKej + MsvjVsl , , 

Mej + Msw ^ ’ 

where Mej is the mass contained in the initial fast moving 
shell at constant velocity I 4 j, and I 4 i the constant velocity 
of the slower moving wind. The mass of the initial slow 
wind swept-up by the shell, Mg^, is given by 

Mgw = ^Rs, (10) 

Fsl 

i.e., the shell mass (Mg = Mej-I-Mg^) will increase linearly 
with radius. 

Assuming I4j=30kms“^ and Mej=8xl0 “^Mq, fur¬ 
ther adopting the median values for the mass-loss rate 
and the gas expansion velocity of a sample of irregular 
and semiregular carbon stars, Msi=3x 10“^ Mq yr“^ and 
t4i=10kms“^ (Schoier & Olofsson 2001) as typical values 
for the slow wind, we hnd that the velocity of the shell 


decrease as illustrated in Fig. 8 . The mass-radius rela¬ 
tion predicted by the model is shown in Fig. 4. In all, this 
simple interacting wind model explains the observed prop¬ 
erties of detached shells, bearing in mind that individual 
objects very likely have different characteristics. 

Under the assumption that the slow stellar wind is 
that of a typical AGB star there is still the need of a very 
brief period of high mass loss also in the interacting wind 
scenario. These eruptive events are most likely linked to 
a He-shell flash (thermal pulse) as originally suggested by 
Olofsson et al. (1990) and Vassiliadis & Wood (1993). 

It is possible that the swept-up material is depleted in 
CO due to it being photodissociated to some degree. This 
would be most pronounced for the older shells. At the 
densities and temperatures prevailing in the shells there 
are no means of effectively producing CO (Mamon et al. 
1988). In this sense the estimated gas shell masses should 
be regarded as lower limits, in particular for the larger 
shells. 

The cooling time scales are ~ 10^ — 10^ yr for gas cool¬ 
ing and as high as 10^ —10®yr for the dust (Burke & 
Hollenbach 1983). Even the gas cooling time scale is rather 
long compared with the shell formation and dynamical 
time scales. However, the time scales will be significantly 
reduced in a clumpy medium. A clumpy medium, with 
gas and dust well mixed, would also explain why the dust 
shells show effects of swept-up material despite the fact 
that the dust mean-free path is rather long for the densi¬ 
ties derived assuming homogeneous shells. 

Finally, we have found evidence that the present-day 
mass-loss rate is decreasing with the size of the shell. The 
most reasonable interpretation of this is that the mass- 
loss rate declines to very low levels after the brief period 
of very high mass loss on a time scale of a few thousand 
years. Thus, there should be material inside the shells. 
There is no indication of this in the interferometer CO line 
maps, suggesting that here the CO molecules have been 
photodissociated. In the shells the densities are higher and 
the CO molecules survive longer. This is supported by the 
detections of ^^CO in some cases. However, it should be 
noted that interferometers are notoriously insensitive to 
weak extended emission and that Lindqvist et al. (1999) 
and Olofsson et al. (2000) only recover about 50% of the 
total flux in their interferometric maps. 

6.3. Photodissociation 

The ^^CO/^^CO-ratio in the shells of V644 Sco and S Set 
are 25 and 20, respectively. For S Set there exist pho- 
tospheric estimates of the ^^C/^^C-ratio. Lambert et al. 
(1986) have estimated this ratio to be 45, significantly 
larger than the value of 14 subsequently obtained by 
Ohnaka & Tsuji (1996) in a different analysis. The rather 
low values for the circumstellar ratio, and in the case of 
S Set a rough agreement with the estimated photospheric 
ratio, suggest that photodissociation of CO molecules is 
low in the shells. If photodissociation would be effective 
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Table 6. Predicted CO line intensities (in Kkms for the 
detached shell around V644 Sco. 


Transition 


TkinjK] 


100 

200 

400 

J = 

1 ^ 0 

6.5 

5.8 

5.1 

J = 

2 ^ 1 

22.6 

22.8 

22.6 

J = 

3^2 

20.7 

22.9 

25.5 

J = 

4^3 

7.8 

12.5 

13.2 

J = 

6^5 

0.58 

1.2 

2.3 

J = 

7^6 

0.15 

0.41 

1.0 

J = 

8^7 

0.04 

0.15 

0.44 


this ratio should be much higher due to the lower ability 
of self-shielding for the rarer isotopomer. The fact 

that the CO molecules appear to survive even out to dis¬ 
tances of « 5x10^”^cm would suggest that the medium is 
clumpy to a large degree, as indeed suggested by observa¬ 
tions (Olofsson et al. 1996; Lindqvist et al. 1999; Olofsson 
et al. 2000). 

6.4. Predictions for future sub-millimetre observations 

To date only three sources have been studied in detail us¬ 
ing interferometry. It is of importance also to map the 
remaining three sources, in particular R Scl. However, 
high-J CO observations are also useful, in particular to 
constrain the kinetic gas temperature in the shell. To il¬ 
lustrate this, predictions for high-frequency CO observa¬ 
tions have been calculated based on the best-fit model for 
V644 Sco using a 15 m telescope. The results are presented 
in Table 6 for three different temperatures in the shell. 
Note that the density has been decreased slightly when 
the temperature was raised in order to maintain a good 
fit to the observed CO lines (J = 1 ^ 0, 2 ^ 1, 3 ^ 2). 

These results show the potential of single-dish tele¬ 
scopes such as CSO, JCMT and the upcoming APEX® to 
further constrain the properties of detached shells. 

7. Conclusions 

We report the detection of a detached molecular shell 
around the carbon star DR Ser. The properties of the 
shell are similar to the young shells found around the car¬ 
bon stars U Cam and V644 Sco. With DR Ser the total 
number of carbon stars with molecular detached shells are 
seven. In fact, there is not much hope of increasing this 
number further since most of the reasonably nearby AGB 

® The Atacama Pathfinder Experiment (APEX), is a collab¬ 
oration between Max Planck Institut fur Radioastronomie (in 
collaboration with Astronomisches Institut Ruhr Universitat 
Bochum), Onsala Space Observatory and the European 
Southern Observatory (ESO) to construct and operate a mod¬ 
ified ALMA prototype antenna as a single dish on the high 
altitude site of Llano Chajnantor. 


stars with mass loss have already been searched for cir- 
cumstellar CO radio line emission, and for the more dis¬ 
tant objects the double-peaked line profiles become less 
pronounced and detached shell sources accordingly more 
difficult to identify. We therefore considered the time ripe 
for a systematic analysis of the existing objects. 

Based on radiative transfer modelling of both observed 
molecular line emission and continuum emission the prop¬ 
erties of the dust and gas in these shells, as well as the 
present-day stellar mass loss have been investigated. It 
turns out that there is a clear trend that both the dust 
and gas shell masses increase with radial distance from 
the central stars. At the same time the velocity by which 
the shells recede from the stars decreases with radial dis¬ 
tance. The most plausible explanation for this behaviour is 
that the shell is sweeping up material from a surrounding 
medium. 

We find that an interacting wind scenario where a brief 
period (~ 10^ yr) of high mass-loss rate (~ 10“® M© yr“^) 
results in a high-velocity shell (;<30kms“^) that expands 
into a previous slow, low-mass-loss-rate AGB wind gives 
an adequate description of the observations. There is also 
evidence that the mass-loss rate, following the short period 
of very intense mass loss, decreases on a time scale of a few 
thousand years to reach a low level (a few 10“® Mq yr“^). 
The most plausible scenario for such a mass-loss-rate mod¬ 
ulation is a He-shell flash (thermal pulse). 

The gas kinetic temperature is the least well con¬ 
strained parameter in the modelling. We suggest that ob¬ 
serving high-J transitions from CO using current and fu¬ 
ture single-dish telescopes should help to constrain this 
parameter. 

Among the seven stars with known detached molecular 
shells, R Scl clearly stands out. The modelling of its ^^CO 
line emission suggests a shell located about a factor of two 
closer to the star compared with observations of scattered 
stellar light. Also the derived mass is significantly higher 
than obtained from dust and ^®CO line modelling. This 
could possibly indicate that R Scl has had a significantly 
different mass-loss-rate history than the other detached 
shell sources, and this warrants further study, in particular 
through interferometric CO observations. 
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